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NEUTRON YIELD CALCULATION FOR GALLIUM-66 FROM DIFFERENT REACTIONS 

B. M. SAIED 1 , A. A. HUSSEIN 2 , 1. N. IBRAHIM 3 , M. N. SARSAM 4 , & R. F. RABEE 5 



The calculated neutron yields are very important in analyzing radiation shielding of spent fuel storage, transport 



and safe handling. The cross sections of ^Zi^p, n)^Ga, f lZn(_p } 2 n) ^Ga, HZn(p } 3n) ^Ga, ¥*Cu (a, 3n) ^Ga, 
f^Zn (a , n -f- p) ^Ga , ^Zi^d, 2n) ^Ga, *2* Cu ( a ; x ) ai Ga and "ao Ell (^ x ) ai Ga are calculated for different 



energies using different sets of programs using Matlab language. The weight average cross section was then used to 
calculate the neutron yields .The empirical formula was then suggested for calculate total neutron yield for each reaction. 
Different type of expressions have been tried to fit the weighted average (W.A) of n-yield data. For the first time new 
polynomial expressions were found relating n-yields with a-energy. 
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The problems of nuclear technology are many: one of them is the development of analytical methods for the 
control and protection from nuclear radiation; another is the production of neutron and isotopic energy sources based on 
radiation emitting of radionuclides. This kind of information must be known to an accuracy of approximately 10% for 

particle energies ranging from reaction thresholds to 10 MeV for target nuclei with z^20 [1] . Medical applications of 

nuclear radiation are of considerable interest to the IAEA. Cyclotrons and accelerators, available in recent years are in 
increasing number in different countries, they are being used for the production of radioisotopes for both diagnostic and 
therapeutic purposes. The physical basis of this production is described through interaction of charged particles, such as 
protons, deuterons and alphas, with matter. These processes have to be well understood in order to produce radioisotopes 
in an efficient and clean manner. In addition to medical radioisotope production, reactions with low energy charged 
particles are of primary importance in nuclear astrophysics there is interest in numerous reaction rates to understand the 
nuclear reactions in the cosmos [2 \ Elemental gallium is not found in nature, but it can obtained by different nuclear 
reactions. Gallium is used in electronics, in microwave circuits and high-speed switching circuits, in infrared circuits in 
blue and violet light-emitting diodes and in diode lasers. The cross section data of g^Zn(p J n^^Ga are obtained from the 

works published by F. Szelecsenyi.et.al [3] in the energy range (25.3- 99.5) MeV, F.Szelecsenyi.et.al [4] in the energy range 
(6.38- 25.72) MeV, B.S.Ishkhanov.et.al [5] in the energy range (4.78- 29.92) MeV, N.A.Demekhina.et.al [6] in the energy 
range (6.4- 21.1) MeV and F.E.Little.et.al [7] in the energy range (6- 35) MeV, these data are plotted, interpolated and 
recalculated in steps of 0.01 MeV using different sets of Matlab programs in order to obtain a weighted average value for 
this reaction. The cross section data of g^Zn^p, 2n)|iGa are obtained from the works published by F. Szelecsenyi.et.al [4] 
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in the energy rang (14.73- 25.71) MeV, A. Henrmannen m in the energy rang (17.94- 29.92) MeV, B.S.Ishkhanov.et.ar J in 
the energy range (13.8- 29.5) MeV and O.V.Ogdankevich.et.al [9] in the energy range (13.6- 21.3) MeV, these data are 
plotted, interpolated and recalculated in steps of 0.01 MeV using different sets of Matlab programs in order to obtain a 
weighted average value for this reaction.The cross section data of I^Znfp, 3n)|^Ga are obtained from the works 

published by F. Szelecsenyi.et.al [10] in the energy range (24.78- 99.26) MeV, T. Stoll et.al [11] , in the energy range (25.7- 
68.9) MeV, A.Hermanne.et.al [12] in the energy range (24.2- 33.9) MeV and T.E.Boothe.et.al [13] in the energy range (25.75- 
29.92) MeV, these data are plotted, interpolated and recalculated in steps of 0.01 MeV using different sets of Matlab 
programs in order to obtain a weighted average value for this reaction. The cross section data of ||C U (a, 3ll) f ^Ga are 

obtained from the works published by V.N. Levkovskij[ 14] in the energy range (28.2- 46) MeV, A. Navin.et.al [15] in the 
energy range (29.1- 37.8) MeV, K. G. Porges [16] in the energy range (29.6- 37) MeV, these data are plotted, interpolated 
and recalculated in steps of 0.01 MeV using different sets of Matlab programs in order to obtain a weighted average value 
for this reaction. The cross section data of f^Zn (a , n + p) Ga are obtained from the works published by V. N. 

Levkovskij [14] in the energy range (18.4- 46) MeV, F. H. Ruddy and B. D. Pate[17] in the energy range (19.6- 36.9) MeV], 
and N.T. Porile [18] in the energy range (18.9- 39.6) MeV, these data are plotted, interpolated and recalculated in steps of 
0.01 MeV using different sets of Matlab programs in order to obtain a weighted average value for this reaction. The cross 
section data of |^Zn{d J 2n)^Ga are obtained from the works published by J. Steyn and B.R. Meyer [19] in the energy 

range (9.7- 15.8) MeV, D.C. Williams and J.W. Irvine Jr [20] in the energy range (9.3- 15.4) MeV and J. L. Gilly.et.al [21] in 
the energy range (8.8- 11.6) MeV, these data are plotted, interpolated and recalculated in steps of 0.01 MeV using different 
sets of Matlab programs in order to obtain a weighted average value for this reaction. The cross section data of 
"S^Cufa^l^Ga are obtained from the works published by F.Szelecsenyi.et.al [22] in the energy range (8.5- 36.3) MeV, 

F.Szelecsenyi.et.al [23] in the energy range (15.9- 58.1) MeV, F.Tarkanyi.el.al [24] in the energy range (7.5- 37.3) MeV, 
M. Sonck.et.al [25] in the energy range (9.9- 42.304) MeV and F.Tarkanyi.et.al [26] in the energy range (8.5- 40) MeV, these 
data are plotted, interpolated and recalculated in steps of 0.01 MeV using different sets of Matlab programs in order to 
obtain a weighted average value for this reaction. The cross section data of "I^Zn^x^^Ga are obtained from the works 

published by F.S.A1-Salenh.et.al [27] in the energy range (3.7- 27.2) MeV, M. S. Uddin et.al [28] in the energy range (4- 39.6) 
MeV, F.Szenlencsennyi.et.al [29] in the energy range (6.42- 95.34) MeV and A. Hermanne [30] in the energy range 
(4.78- 29.29) MeV, these data are plotted, interpolated and recalculated in steps of 0.01 MeV using different sets of Matlab 
programs in order to obtain a weighted average value for this reaction. Knowledge of total neutron yield of heavy ions 
from thick target is important for estimating radiation damage by lost beam ions in high intensity accelerators. The cross 
sections used (o L ) are the weighted average cross sections and their corresponding errors (Atfj) calculated according to the 
following expressions [31] : 

S.D= , 1 , (2) 
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The yield of neutrons per each incident particle Y n for an ideal, thin and uniform target and non energetic beam of 
energy E is given by [32] 

m>£^ -..(3) 

dE 

For each energy interval (E b -E th ), the ratio of the cross section values by the energy loss — {h) calculated using 

the stopping power SRIM gives the yield for that interval, E th is the reaction threshold energy and E b is the average 
energy between two successive values. 

RESULTS 

• The reaction 30 En (ft n )31 Ga 

The cross section data measured by [3-7], have been plotted, interpolated, and recalculated in steps of 0.01 MeV 
from energy 5 MeV up to 50 MeV in order to obtain the adopted cross section of this reaction. The results of adopted cross 
section and the actual data are shown in figure 1 the weighted average cross section can be expressed by the following 
polynomial: 

f(x)=439.1*exp(-((x-16J 128.9*exp (-((x-43.35)/11.41) A 2) 

The obtained neutron-yield using SRIM codes in the chosen energy range are shown in figure 9. 

• The reaction 30 ^ ■ /31 

The cross section data of these reactions measured by [4, 5, 8, and 9], have been plotted, interpolated, and 
recalculated in steps of 0.01 MeV from energy 13 MeV up to 30 MeV in order to obtain the adopted cross section of this 
reaction. The results of adopted cross section and the actual data are shown in figure 2 the weighted average cross section 
can be expressed by the following polynomial: 

f(x) = 618*exp(-((x-29.3)/6.524) A 2) +446*exp(-((x-37.3)/9.249) A 2) 

The obtained neutron-yield using SRIM code in the chosen energy range are shown in figure 10 

• The reaction 30 ^ y31 

The cross section data measured by [16-19], have been plotted, interpolated, and recalculated in steps of 0.01 
MeV from energy 24.5 MeV up to 70 MeV in order to obtain the adopted cross section of this reaction. The results of 
adopted cross section and the actual data are shown in figure 3 the weighted average cross section can be expressed by the 
following polynomial: 

F(x) = -0.4995*x A 4 + 19.47*x A 3 -264.3*x A 2+1497*x-2948 

The obtained neutron-yield using SRIM code in the chosen energy range are shown in figure 1 1 

. Thereaction2l Cu t a ' 3ll )ii Ga 

The cross section data measured by [14-16], have been plotted, interpolated, and recalculated in steps of 0.01 
MeV from energy 27.5 MeV up to 46 MeV in order to obtain the adopted cross section of this reaction. The results of 
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adopted cross section and the actual data are shown in figure 4 the weighted average cross section can be expressed by the 
following polynomial: 

F(x) = 2.303e+4.1*exp(-((x-15.25)/3.8) A 2)+55.66*exp(-((x-9.425)/2.287) A 2)- 2.294e+4.05*exp (-((x-3.793) 
A 2) +34.96*exp (-((x-36.04)/27.58) *2) 

The obtained neutron-yield using SRIM code in the chosen energy range are shown in figure 12 

• The reaction f ^ Zn (a , n + p) f f Ga 

The cross section data measured by [14, 17, 18], have been plotted, interpolated, and recalculated in steps of 0.01 
MeV from energy 20 MeV up to 40 MeV in order to obtain the adopted cross section of this reaction. The results of 
adopted cross section and the actual data are shown in figure 5 the weighted average cross section can be expressed by the 
following polynomial: 

f(x) = 1.916e+4*exp(-((x-14.38)/4.328) A 2)- 1.871e+4*exp(-((x-14.4)/4.303) A 2) +57.26*exp(-((x-26.87)/15.28) A 2) 
The obtained neutron-yield using SRIM code in the chosen energy range are shown in figure (13). 

• The reaction 3 ° 31 - 

The cross section data measured by [19-21], have been plotted, interpolated, and recalculated in steps of 0.01 
MeV from energy 9 MeV up to 15 MeV in order to obtain the adopted cross section of this reaction. The results of adopted 
cross section and the actual data are shown in figure 6 the weighted average cross section can be expressed by the 
following polynomial: 

f(x) = 244.4*exp(-((x-19.9)/2.685) A 2) + 418.7*exp(-((x-23.69)/4.422) A 2) + 183.0*exp (-((x-16.6)/1.853) A 2). 
The obtained neutron-yield using SRIM code in the chosen energy range is shown in figure (14). 

• The reaction "°9 Cu (^ x ) 3i Ga 

The cross section data measured by [22-26], have been plotted, interpolated, and recalculated in steps of 0.01 
MeV from energy 8 MeV up to 58 MeV in order to obtain the adopted cross section of this reaction. The results of adopted 
cross section and the actual data are shown in figure 7 the weighted average cross section can be expressed by the 
following polynomial: 

F(x)= 127.4*exp(-((x-38.26)/1.687) A 2)+153.5*exp(-((x-36.38)/4.824) A 2) 

The obtained neutron-yield using SRIM code in the chosen energy range are shown in figure 15 

• The reaction 30 w J31 

The cross section data measured by [27-30], have been plotted, interpolated, and recalculated in steps of 0.01 
MeV from energy 5 MeV up to 30 MeV in order to obtain the adopted cross section of this reaction. The results of adopted 
cross section and the actual data are shown in figure 8 the weighted average cross section can be expressed by the 
following polynomial: 
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f(x) = 7.791*exp(-((x-33.29)/0.4236) A 2) + 4.05e+15*exp(-((x+655.9)/125.3) A 2)- 625.7*exp(-((x-22.72)/8.363) A 2) 
The obtained neutron-yield using SRIM codes in the chosen energy range are shown in figure (16). 
CONCLUSIONS 

The weighted average cross sections values indicate clearly the necessity to adopt such calculations which are 
very important for some problems such as neutron dosimetry, fuel burn-up and determination of isotope production. 
A complete dependency on individual or even collective experimental results are not recommended due to experimental 
deviations and errors and the impossibility of measurement for all the detailed interval of energy .This may be one of the 
important reasons to explain the observed deviations in some energy intervals of the international libraries. Thus using the 
weighted average cross sections one can obtain an accurate, complete and energy detailed cross section values and can 
implement the essential condensation calculation in terms of any energy interval which are essential for all kinds of reactor 
calculations. By using these adopted cross sections values together with stopping power from (SRIM) software the neutron 
yields empirical formulas were suggested for each reaction, where these polynomial expressions have been established and 
found to be most adequately fits the data. 




Figure 1: Cross sections of 66 Zn (p, n) 66 Ga Figure 2: Cross sections of 67 Zn (p, 2n) 66 Ga 




Figure 5: Cross Sections of 64 Zn (A, N+P) ^ Figure 6: Cross Sections of 66 Zn (D, 2n) 66 Ga 



www.tjprc.org 



editor@tjprc.org 



32 



B. M. Saied, A. A. Hussein, I. N. Ibrahim, M. N. Sarsam & R. F. Rabee 




Figure 7: Cross sections of °Cu (a, x) 66 Ga Figure 8: Cross sections of °Zn (p, x) 66 Ga 




Energy (MeV) 

Figure 9: Total Neutron yield of 66Zn (p, n) 66 Ga Reaction 
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Figure 10: Total Neutron yield of 67Zn (p, 2n) 66 Ga Reaction 




Figure 11: Total Neutron Yield of bbLn (P, 3n) 66 Ga Reaction 
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Figure 12: Total Neutron Yield of 65Cu (A, 3n) 66 Ga Reaction 
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Figure 13: Total Neutron Yield of 64Zn (A, N+P) b0Ua Reaction 
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Figure 14: Total Neutron Yield of 66Zn (D, 2n) 66 Ga Reaction 
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Figure 15: Total Neutron Yield of Nat cu (A, X) 66 Ga Reaction 
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